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1 Introduction. 

Recent years have witnessed some progress in theoretical description of heavy quark systems 
near the threshold which will be experimentally studied at future colliders ]]]. It is mostly a 
technical progress within general approach of effective theories 0. The computational difficul- 
ties of high order PT calculations have led to the necessity of introducing some approximations 
realized in the framework of the effective theories. Processes with a low momentum transfer for 
the particles containing one heavy quark are well approximated by HQET ||. This approach 
was quite successful in describing transitions between heavy particles containing b- and c-quarks 
0| . For the systems of two heavy particles near the threshold the effective theory is NRQCD 
. Recent result obtained in its framework is the completion of NNLO analysis of heavy quark 
production e.g. ||. In my talk I briefly present the results of computing perturbative correc- 
tions to physical quantities relevant for description of baryons with one heavy quark and of the 
heavy quark pair production near the threshold. 

2 Baryons with one heavy quark. 

Properties of the spectrum of heavy baryons can be obtained from the analysis of the correlator 
of two baryonic currents. Calculations with massive quarks are basically done in the leading 
order only. The NLO corrections in the massless case have been known since long ago 0. The 
generic baryonic current is a local three-quark operator j = e abc (uJCd c )T^/ a where \1/ is a heavy 
quark with the mass m, u, d are massless quarks, C is the charge conjugation matrix, T is a 
Dirac matrix. e abc is the totally antisymmetric tensor, a, b, c are indices for the SU(3) color 
group. Both functions U qtTn (q 2 ) of the correlator of two baryonic currents (for simplicity we 
take T = 1) 

if(T 3 (x)m)z iqX dx = lu q v Yl q (q 2 )+mYl m (q 2 ) (1) 

are known up to NLO of PT (three- loop diagrams). I shall present results only for the function 
U m (q 2 ) [§]. For the spectral density p(s) in the dispersion relation 

IT (a 2 ) - — H ^^ (2) 

we write 

P(s) = s 2 (p (s)(l + * In 4) + -Pi(s)) . (3) 

Here \i is the scale parameter, a s = a s (/i), m is a pole mass of heavy quark 0. The leading 
order term reads 

Po ( s ) = l + g z - g z 2 - z 3 + 6z(l + z) In z (4) 



with z = m 2 /s. In the MS-scheme the NLO term is given by 

/ x r, 665 665 2 n o .„. 

Pi(«) = 9 + — *-— z 2 -^ 3 ( 5 ) 
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- (y + 42z - 42^ 2 - ^ 3 ) ln(l - s) + (2 

+4 (3 + 3^ - 3^ 2 - -r 5 ) ln(l - z) lnz + 12z (2 + 3z + -z 2 ) (- In 2 z - ((2) 
+4 (- + Viz + 3z 2 - -z 3 ) U 2 (z) + 24z(l + z) ^{z) - C(3) - -Li 2 (z) lnz) 
where Li n (-?) are polylogarithms and ^(n) is Riemann's zeta function. The diagrams have 



generalized sunset topology which has been recently studied in some detail |LC| . Such a topology 
of diagrams is rather frequent in phenomenological applications JD]]. Since the anomalous 
dimension of the baryonic current is known to two-loop order [|12|] the result eq. (j|) completes 
the NLO analysis. 

Two limiting cases of interest are the near-threshold and high-energy asymptotics. With 
eq. (H) both limits can be taken explicitly. Effective theories can be viewed as special devices 
for calculations in such limiting cases. 

In the high energy (small mass) limit z — > the correction reads 

Pl ( s ) = 9 + 83z-4n 2 z + 2 ln(Y) + 50,2 ln(z) + 12 z\n(z) 2 - 24 z( (3) + 0{z 2 ) . (6) 
This leads to the small mass expansion of the spectral density in the form 



mp(s) = m m (fi)p massless (s) = m m (p)s 2 1 + — 2 In — + — (7) 



where p m assiess(s) is the result of calculating the correlator in the effective theory of massless 
quarks. The relation between the pole mass m and the MS mass fnj^(p) reads 



7T V 77T 3 



m = mMs(/x)(l + ^(ln^ + -)) . (8) 



The massless effective theory cannot reproduce the mass singularities - terms like z ln(z) in 
eq. (||). Within the massless effective theory such contributions can be parametrized with 
condensates of local operators. The first m 2 correction in eq. @ can be found if the pertur- 
bative value of the heavy quark condensate (\I/\I/) taken from the full theory is added. The 
vacuum expectation value (\I/\I/) of the composite operator (\E r \l / ) should be understood within 
a mass independent renormalization scheme such as the MS-scheme. This value (perturbative, 
(ifrty ~ m 3 ln(/i 2 /m 2 )) cannot be computed within the effective theory of massless quarks. It 
represents the proper matching between the perturbative expressions for the correlators of the 
full (massive) and effective (massless) theories. This matching procedure allows one to restore 



higher order terms of the mass expansion in the full theory from the effective massless theory 
with the mass term treated as a perturbation. 

In the near-threshold limit E — * with s = (m + E) 2 one explicitly obtains 

M",S) = — (l + -ll„^ + -i( y + _ + 41n-)) + 0(^). (9) 

In this region the appropriate device to compute the limit of the correlator is HQET. Writing 

mpth{m, E) = C(m/p, a s ) 2 p HQET (E, p) (10) 

we obtain the known result for p H q ET (E, p) |L3| 

/Tn N 16E 5 { a s /182 4vr 2 ,, u \\ _ ., , . 

P„ Q et(E,») = — (l + f — + — + 41n^ + 0(E") (11) 



with matching coefficient C{mj 'p, a a ) [|T4 
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The matching procedure allows one to restore the near-threshold limit of the full QCD correlator 
([5J) starting from the result obtained in a simpler HQET valid near the threshold only. Higher 
order terms of E/m expansion can be easily obtained from the explicit result eq. fl5|). The NLO 
correction in low energy expansion reads 

A , ^ 88E 6 f a,/, ( p 2 \ 376 4vr 2 140, f m\\\ 

Ap.,KB) = -^{l + ^(l.>^) + - + - + -ln(-)j|. (13) 

To obtain this result starting from HQET is a more difficult task requiring the analysis of 
contributions of higher dimension operators. 

Note that the simple interpolation between the massless and HQET limits gives a rather 
good approximation of the ratio pi(s)/po(s) for all m 2 < s. For the moments of the spectral 
density 

M n = r ^ = m e ~ 2n M n , M n = Mf (l + «• (In 4 + *»)) (14) 

one finds 

Ml 0) = — ,„, 12 ,„, -, S n = A n + -C(2) = A n + — . (15) 

n(n - l) 2 {n - 2) 2 (n - 3)' 3 SV ' 9 V ; 

The coefficients A n are rational numbers, A 4 = 3, v4 5 = 13/2, v4 6 = 17/2, A 7 = 535/54 ||. The 
expression for <5 n with arbitrary n is long. 



3 Two heavy quarks near threshold. 

With two heavy quarks involved in a physical process it is becoming more difficult to perform PT 
calculations. Still many results in finite order PT for vector current correlators are known both 



analytically and numerically ||15| . For applications near the threshold where v = Jl — Am 2 /s <C 
1 (m t is a generic notation for the heavy quark mass) the ordinary perturbation theory (with 
free quarks as the lowest order approximation) breaks down. The ratio a s /v is not small and 
all terms of the order (a s /v) n should be summed. The motion of the heavy quark-antiquark 
pair near the production threshold is nonrelativistic to high accuracy that justifies the use 
of nonrelativistic quantum mechanics for describing such a system ||. Being much simpler 
than the comprehensive relativistic treatment of the bound state problem with Bethe-Salpeter 
amplitude |16[], this approach allows one to take into account Coulomb interaction exactly 
|17| , |i8| . However, for this approximation to work the full QCD should first be mapped onto 
NRQCD which is done in PT. For the analysis of heavy quark production by the vector current 
the basic quantity is the correlator 

U(E)=zJ(Tj em (x)j em me iqx dx } q 2 = (2m t + E) 2 . (16) 

Near the threshold (for small energy E) NRQCD is used. Up to the NNLO accuracy in NRQCD 
one has (cf. eq. fllOl) ) 

11(E) = ^C h (a s )C (E/m t )G(E; 0, 0) . (17) 

Here Ch{a s ) is the high energy coefficient (analog of eq. fll2]) ) which has been known in the NLO 
since long ago |]1|. G(E; 0, 0) is the nonrelativistic Green's function (an object of the effective 
theory as Phqet(E, fi) in eq. ([TT])), E = y/s — 2m t . Up to NNLO of NRQCD the contribution 



of higher dimension operators can be written as a total factor Co(E/m t ) = 1 — AE/3m t [p0 |. 
High energy coefficient Ch(a s ) is given by the expression 

C 4 (a) = l-4^ + C^) 2 (-^ln« + c 2 ). (18) 

The most difficult part of eq. ([H|) to obtain is the coefficient c 2 [|21 
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c 2 = + 7T 2 n 2 In 2 C(3) + —Tit, 
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rif is the number of light quarks. An explicit dependence of high and low energy quantities 

on the factorization scale /i/ is a general feature of effective theories which are valid only for 

a given region of energy. A physical quantity, which is given by a proper combination of 

results obtained in different energy regions within respective approximations, is factorization 

scale independent e.g. [22] . The basic dynamical quantity in the analysis of the near-threshold 



effects is the nonrelativistic Green's function G(E) = (H — E) 1 where 



2 

H=^-+V(r) (19) 

m t 



is the nonrelativistic Hamiltonian. A part of Hamiltonian ( |T9"D which is difficult to find is 
the heavy quark static potential V pot (r) entering into the potential V(r). The static potential 
Vp t{r) is computed in perturbation theory 

VWr) = "CVy (l + ^ (<?i Mjv) + Co) + (^) 2 (Cl InV) + C\ ln(H + C 2 )) . (20) 

Nontrivial coefficients of this expansion are Cq and Cq. The numerical value for Cq was obtain 
long ago |23[, while Cq has been recently computed |2"4]1 . In order to exactly take into account 



the Coulomb effects near the threshold the Hamiltonian ([191) is represented in the form [25| 



2 

H = H C + AH, H C = — -C F — (21) 

m t r 



with 
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r. 

For QCD one has Ca = 3, Cf = 4/3. Constructing the Green's function is straightforward and 
can be done analytically within perturbation theory in AH near the Coulomb Green's function 



Gc(E) |26|, £3, gg] or numerically by solving the Schrodinger equation with Hamiltonian ( |19|) 
(for complex values of E that can be used to describe the production of particles with nonzero 
width) |29], |30], |3l], |32], ^3| [54j. I'll discuss the analytical solution only. Some terms in eq. ( |22"D 
can be handled efficiently. The shift of the parameters of the Coulomb Green's function E — > 
E + E 2 /4m t and a s — > a s (l + 3E/2m t ) accounts for the relativistic H^ and anticommutator 
corrections in eq. (^) . In this respect the modified Coulomb approximation can be used as the 
leading order approximation. The Hamiltonian can be written in the form 

4tt fC F C A \ 2vr 3") 



H = H w + AV M + a ,V a S { r ) , * = -£(£ + %)=*£ (23 ) 



where H LO is the (modified) Coulomb approximation, AV pot is the perturbation theory correc- 
tion to V po t(r). The <5(r)-part is a separable potential and can be taken into account exactly 
The solution reads 



g(^0,0)= 1 ^ G : / ( ^°;° ) nm , G ir {E) = {H LO + AV PT -E)-\ (24) 

1 + a s VoG ir (E; 0, 0) 

Therefore one is left only with the necessity to construct corrections due to AV^ oi -part which 
modify Gi r (E). An efficient calculational framework for perturbation theory in AV pot near the 
Coulomb Green's function Gc{E) has been developed p6| . The key point is the use of a partial 
wave decomposition for the Green's function 

oo 

G(E; x, y) = ]T(2/ + l){xy) l P l ^y/xy)G l {E ] x, y) (25) 

1=0 
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where Pi(z) is a Legendre polynomial. The Coulomb partial waves Gf{x, y, k) are known 

Ll + \2kx)Ll + \2ky)m\ 



■''"■■■ ■■ - - mtk (2k) 2l e- k ^ x+ ^ J2 



Gf(E;x,y) 



2ir 



'„ (m + J + 1 - i/)(m + 2/ + 1)! 



(26) 



where fc 2 = —m t E, v = X/k, X = a s Cpm t j2, L^(z) is a Laguerre polynomial. The result of 
the evaluation of the NLO correction in S-wave is 

A.GoiE) = -£^—± £ F(m) 2 (m + 1) (L^k) 4- ^(mH- 2)) - 2 £ £ F(m) 

Z7r Z7r \m=0 m=l n=0 

r? 4- 1 °° 1 \ 

XjP ( n )^^_ + 2 ^ F(m) (Li(A;) - 2 lE - ^(m + 1)) - 7^1 (*0 + 7^1 (*0 2 
where *n(aO = <i n lnr(a;)/<ia; n , /?o is the first coefficient of the /5-function. Here 



Li(fc) =ln 



2A- 



F(m) 



(m + 1) (m 4- 1 — u) 



The NLO correction to the I = 1 partial wave was found in |36 
a s j3o Xm t k 2 ' v 



A l G 1 (E) 



2tt 18tt 



J2 F(m) 2 (m 4- l)(m 4- 2)(m 4- 3) (L X (A;) + *i(m 4- 4)) 



Vm=0 



00 m— 1 



-2 X! E ^)^W (n + 1)(n + 2)(n + 3) + 2 f; F(m) ( 2J (m) + (m + l)(m + 2)£ 1 (ib) 



m=l n=0 



m — n 



m=0 



+ (1 4- z/)(Ji(m) + (m 4- l)Li(fc)) + ^-^ (J 2 (m) + 2L 1 (A;)) ) 4- /(fc) 



where 



F(m 



v(v 2 - r 



(m 4- 2 - i/)(m 4- l)(m 4- 2)(m + 3) 
Jo(m) = -2*i (m + 1) - 4 7B + 3, 

J^m) = (m 4- l)(-*i(m + 2) - 2 7s + 2), 

~ . , (m+l)(m + 2) / _ , , n 3^ 

^(m) = ^ ^ '- f -*!(m + 3) - 2 7B + - , . 



(7fl ~ I) 2 ^ 2 ,„ o n , 1 - 9 7 j; 4- 67J 4- 7T 2 2 l-37g 3 1-tj; 4 
4 — 37r )v H z/ H 1/ H z/ 
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+ I 7^ - 1 — 3i/ • 



9 " 127 "- 2 + ?, 3 + -A Ll (k) +(- l - + l" 2 ) Lr{kf 



I(k) 



4 2 4/ v/ V22 

The Green's function at the origin can be written in the form with single poles only 

1 ; ^ E 0m -E Wo E'-E 



(27) 



where Vw.^O) is the wave function at the origin. Up to NNLO one writes 



Eom 
2 



E^ m (l + A 1 E 0m + A 2 E l 



'OmJ 



IVw(o)| 



|^o C m (0)| 2 (H-A^ 2 m + A 2 ^ 2 m 



(28) 



where 



E, 



c 

(hit 



X 2 



m t (m + 1) 



IV>LW 



A :3 



7r(m + l) 3 



The analytical expression for the NLO corrections to the bound state parameters has the form 

a s p 



A 1 E I 



Orn 



TX 



(L x {m) + *i(m + 2] 



a s p 



AiV'om = : ^p(3^i(m) + ^ 1 (m + 2)-2(m + l)^ 2 (m+l)-l-27 £ + 
where Li(m) = Li(\/(m + 1)). The expressions of the NNLO corrections to the energy levels 



and wave functions at the origin are rather long |27], |37], |38fl . As an example I give the correction 
A 2 Vom due to the second iteration of the 0{a s ) correction to the static potential P2 



,W./,2 



a 



.V 



-A"'' Vi„. = i ^J {3(Ci + (Hm) + Mm + 1))C\f 



+2Cl ft C. + i)( m +i) ^ + / £ + + i. » + 1 



2 (n — m) (m + 1) 



^ 



V { m + 1} ' fc 1 + f (Z(m) + *i(n + 2) - -^- ) ( <\ 
^ (n-m) 3 V ° V V V ; V ; 2n-m' 



n=m+l 



+ 2C{ (C^ + (Z(m) + *!(m + 2)) C?) -- + £ 



5 ^ n + 1 



n=0 



(n — m) 



;U(m,n) 



^^^-"jj^^li-s 



I ^ n + 1 



(n — m) 2 



1 ^ n + 1 



+^E 



2 „=o(«- ™) 3 



„, v 1 ^ (m + i) 



n=m+l 



(n — m) 3 (n + 1) 



~ m+1 

„=™+i (n - m) 2 

U(m,n) 



(Z + l)(m + l) 



(/ + l)(m + l) 



, yyV (Z + l)(n+l) 

~[ ^ \(n — m) 2 (/ — m) 2 (n — m) 2 (l — m)(n — 1) (I — m) 2 {n — m)(n — I) 



~ ^V (* + l)(m+l) (l + l)(m + l) 2 

n^m+i z^o \ {n- m) 2 (l - m) 2 (n - m) 2 (l - m)(n - l)(n + 1) 



(i + l)(m + l) \ + ~ ^ / (m + 1) 2 

(/-m) 2 (n-m)(n-/)/ ^J^ \(n - m) 2 (l - m) 2 



(7 + l)(m+l) 2 (m + 1) 

"I" / \0/T \/ 7\ / . -. \ 



(n — m) 2 (l — m)(n — l)(n + 1) (/ — m) 2 (n — m)(n — /) 
where L(m) = L(\/(m + 1)) and 

n+1 (n + 1) 2 

U(m,n) = 3+ — -^— r -2- 



m + n + 2 (n — m) (n + m + 2) 

4 Physical results. 

There are basically two application areas of the near-threshold calculations within NRQCD: the 
analysis of T-resonances in the bb system and the ti production cross section near the threshold. 

4.1 T-resonances. 

For the bb system moments of the spectral density of different kinds can be studied (cf. eq. (|14D). 

For instance, 

p(s)ds 



/w s 11 

b 

are meaningful in the near-threshold Coulomb PT calculations [p9 ]. In the NNLO we obtain 
the numerical value for the 6-quark pole mass p] 



m b = 4.80 ± 0.06 GeV. (29) 

With this result we obtain the value of the matrix element \V c b\ from the analysis of the 5-meson 
semileptonic width [flTj 

I*, - 0.0.3 (^=i^f (^f (l - 0M^_1^) (1 ± ^ (30) 

where A npt ~ 0.02 is the uncertainty of nonperturbative contributions. We use the analysis 
in the pole mass scheme. Other approaches based on redefinition of the mass near threshold 
have been suggested. The problem of infrared instability of the pole mass and its relation to 
MS-mass has been discussed [^7[ [R| . 

4.2 ti cross section near threshold. 

The it-pair near the production threshold is just a system that satisfies the requirement of being 
nonrelativistic. Therefore the description of tt-system near the production threshold y/s w 2mt 
(\fs is a total energy of the pair, m t is the top quark mass) is quite precise within NRQCD. The 



top quark is very heavy m t = 175 GeV [41] and there is an energy region of about 8 — 10 GeV 



near the threshold where the nonrelativistic approximation for the kinematics is very precise. 

9 



Relativistic effects are small and can be taken into account perturbatively. The strong coupling 
constant at the high energy scale is small a s (m t ) « 0.1 that makes the mapping of QCD onto 
the low energy effective theory (NRQCD), which is perturbative in a s (m t ), numerically precise. 
The top quark decay width is large, T t = 1.43 GeV; the infrared (small momenta) region is 
suppressed and PT calculations for the cross section near the threshold are reliable even point- 
wise in energy E2|. One can study the ti system near the threshold in the processes e + e~ —>■ ti 



and 77 — > tt [43 



e + e — ► tt: the production vertex is local, the basic observable is a production cross 
section which is saturated by S-wave (for the vector current), NNLO analysis is available. 

77 — > ti: the production vertex is nonlocal (T-product of two electromagnetic currents), 
both S- and P-waves can be studied for different helicity photons, the number of observ- 
ables is larger (cross sections as, ap, S-P interference). Because of nonlocality of the 
production vertex the high energy coefficient (necessary for mapping the QCD quantities 



to NRQCD ones) is more difficult to obtain. Some calculations were done in NLO [44 , [45 



and the full analysis of the cross section is available in NLO of NRQCD only |36|. The 



low energy part of the process can be studied in NNLO without a strict normalization to 
full QCD (see [46| for relativistic corrections). 



The top quark width T t plays a crucial role in the calculation of the tt production cross section 
near the threshold. It is taken into account by the shift E — > E + iY t ||42|| . One has 

p{E')dE' ^ f p{E')dE' 

{E'-Ef + Yl 
Because the point E + iYt lies sufficiently far from the positive semiaxis in the complex energy 



^^nM^jJ^-r,/*™*. (31) 



plane the cross section eq. (^1|) is calculable point-wise in energy. The hadronic cross section 
a(E) was obtained by many authors (as a review see, ||). The normalized cross sections R V (E) 
are plotted in Fig.l | |32j j . The curves have characteristic points which are usually considered 
as basic observables. They are: E p - the position of the peak in the cross section and H p - 
the height of the peak. The convergence for E p and H p in consecutive orders of PT near the 
Coulomb solution is not fast in the MS-scheme [|, [38], f37| . For typical values m t = 175 GeV, 
T t = 1.43 GeV, a s (M z ) = 0.118 one finds pi 



E p = Eo(l + 0.58 + 0.38 + ...) 

H p = H (l - 0.15 + 0.12 + ...) . (32) 

The way of dealing with the PT expansion of the static potential is an important issue in 
getting stable results within NRQCD [|4"8" |. The convergence in the MS-scheme is not fast which 



reflects the physical situation that the observables E p and H p are sensitive to different scales. 
The finite-order perturbation theory expansion of the static potential given in eq. (|2T)| ) cannot 
handle several distinct scales with the same accuracy. 

10 




Fig. 1. The normalized cross section R V (E) in the leading order (solid lines), NLO 

(bold dotted lines) and NNLO (bold solid lines) for m t = 175 GeV, T t = 1.43 GeV, 

a s (M z ) = 0.118 and /i s = 50 GeV, 75 GeV and 100 GeV. The dotted line corresponds 

to the result in Born approximation. 

5 Conclusion. 

Some calculations of PT corrections to systems with heavy particles are presented. The use 
of effective theories for computing high order corrections is crucial in obtaining precise results 
near the threshold. The calculation are rather tedious but still can be completed analytically 
in some cases. 

Acknowledgments 

This work is partially supported by Volkswagen Foundation under contract No. 1/73611 and 
Russian Fund for Basic Research under contract 99-01-00091. 



References 

[1] E. Accomando et al, Phys. Rep. 299 (1998) 1. 

[2] S. Weinberg, Physica A96 (1979) 327; J. Gasser, H. Leutwyler, Nucl. Phys. B250 (1985) 
465. 



11 



[3] M.A. Shifman, M.B. Voloshin, Sov. J. Nucl. Phys. 45 (1987) 292; H.D. Politzer, M.B. 
Wise, Phys. Lett. B206 (1988) 681, Phys. Lett. B208 (1988) 504; E. Eichten and B. Hill, 
Phys. Lett. B234 (1990) 511; H. Georgi, Nucl. Phys. B363 (1991) 301; M. Neubert, Phys. 
Rep. 245 (1994) 259. 

[4] N. Isgur, M.B. Wise, Phys. Lett. B232 (1989) 113, Phys. Lett. B237 (1990) 527. 

[5] W.E. Caswell and G.E. Lepage, Phys. Lett. B167 (1986) 437; G.E. Lepage et al, Phys.Rev. 
D46 (1992) 4052; G.T. Bodwin, E. Braaten and G.P Lepage, Phys.Rev. D51 (1995) 1125. 

[6] A.H. Hoang et al, Eur. Phys. J.direct C3 (2000) 1. 

[7] A.A. Ovchinnikov, A. A. Pivovarov, L.R. Surguladze, Sov. J. Nucl. Phys. 48 (1988) 358, 
Int. J. Mod. Phys. A6 (1991) 2025. 

[8] S. Groote, J.G. Kdrner and A.A. Pivovarov, Phys. Rev. D61 (2000) 071501. 

[9] R. Tarrach, Nucl. Phys. B183 (1981) 384. 

[10] F.A. Berends, A.I. Davydychev, N.I. Ussyukina, Phys. Lett. B426 (1998) 95; S. Groote, 
J.G. Kdrner and A.A. Pivovarov, Phys. Rev. D60 (1999) 061701, Eur. Phys. J. Cll (1999) 
279, Nucl. Phys. B542 (1999) 515, Phys. Lett. B443 (1998) 269; S. Groote and A.A. Pivo- 
varov, Nucl. Phys. B580 (2000) 459; A.I. Davydychev and V.A. Smirnov, Nucl. Phys. 
B554 (1999) 391; N.E. Ligterink, Phys. Rev. D61 (2000) 105010. 

[11] J.O. Andersen, E. Braaten, M. Strickland, Phys. Rev. D62 (2000) 045004; S. Narison and 



A.A. Pivovarov, Phys. Lett. B327 (1994) 341; T. Sakai, K. Shimizu and K. Yazaki, [nucl 
th/9912063j ; S.A. Larin et al, Sov. J. Nucl. Phys. 44 (1986) 690; J.M. Chung, B.K. Chung, 
Phys. Rev. D60 (1999) 105001; K. Chetyrkin, S. Narison, Phys. Lett. B485 (2000) 145; 
H.Y. Jin, J.G. Kdrner, |hep-ph/0003202| . 

[12] A.A. Pivovarov, L.R. Surguladze, Yad. Fiz. 48 (1988) 1856; Nucl. Phys. B360 (1991) 97. 

[13] S. Groote, J.G. Kdrner and O.I. Yakovlev, Phys. Rev. D55 (1997) 3016. 

[14] A.G. Grozin and O.I. Yakovlev, Phys. Lett. 285B (1992) 254. 

[15] K.G. Chetyrkin, J.H. Kuhn, A. Kwiatkowski, Phys. Rept. 277 (1996) 189; K.G. Chetyrkin, 
J.H. Kuhn, M. Steinhauser, Phys. Lett. B371 (1996) 93, Nucl. Phys. B482 (1996) 213. 

[16] E.E. Salpeter, H.A. Bethe, Phys. Rev. 84 (1951) 1232. 

[17] M.B. Voloshin, Nucl. Phys. B154 (1979) 365; Yad. Fiz. 29 (1979) 1368, Sov. J. Nucl. 
Phys. 36 (1982) 143; H. Leutwyler, Phys. Lett. B98 (1981) 447. 

12 



[18] W. Kwong, Phys. Rev. D43 (1991) 1488; M.J. Strassler and M.E. Peskin, Phys. Rev. D43 
(1991) 1500; M. Jezabek, J.H. Kuhn and T. Teubner, Z. Phys. C56 (1992) 653; Y. Sumino 
et ai, Phys. Rev. D47 (1993) 56. 

[19] R. Karplus and A. Klein, Phys. Rev. 87 (1952) 848; G. Kallen and A. Sarby, K. Dan. 
Vidensk. Selsk. Mat.-Fis. Medd. 29 (1955) , N17, 1; R. Barbieri et ai, Phys. Lett. B57 
(1975) 535; I. Harris and L.M. Brown, Phys. Rev. 105 (1957) 1656; R. Barbieri et ai, 
Nucl. Phys. B154 (1979) 535. 

[20] A.H. Hoang, Phys.Rev. D56 (1997) 7276. 

[21] A. Czarnecky and K. Melnikov, Phys. Rev. Lett. 80 (1998) 2531; M. Beneke, A. Signer 
and V.A. Smirnov, Phys. Rev. Lett. 80 (1998) 2535. 

[22] A.A. Pivovarov, Phys. Lett. B236 (1990) 214, Phys. Lett. B263 (1991) 282. 

[23] W. Fisher, Nucl. Phys. B129 (1977) 157; A. Billoire, Phys. Lett. B92 (1980) 343. 

[24] M. Peter, Phys. Rev. Lett. 78 (1997) 602; Nucl. Phys B501 (1997) 471; Y. Schroder, Phys. 
Lett. B447 (1999) 321. 

[25] A.I. Achieser and V.B. Berestezki, Quantum electrodynamics, Moscow 1959; L.D. Landau 
and E.M. Lifshitz, Relativistic Quantum Theory, Part 1 (Pergamon, Oxford, 1974) . 

[26] J.H. Kuhn, A.A. Penin and A.A. Pivovarov, Nucl. Phys. B534 (1998) 356. 

[27] A.A. Penin and A.A. Pivovarov, Phys. Lett. B435 (1998) 413; Phys. Lett. B443 (1998) 
264; Nucl. Phys. B549 (1999) 217. 

[28] M. Beneke, A. Singer and V.A. Smirnov, Phys. Lett. B454 (1999) 137. 

[29] A.H. Hoang and T. Teubner, Phys.Rev. D58 (1998) 114023. 

[30] K. Melnikov and A. Yelkhovsky, Nucl. Phys. B528 (1998) 59. 

[31] T. Nagano, A. Ota and Y. Sumino, Phys. Rev. D60 (1999) 114014. 

[32] A.A. Penin and A.A. Pivovarov, MZ-TH-98-61, Dec 1998. 41pp. [|iep-ph/9904278 |. 

[33] O. Yakovlev, Phys. Lett. B457 (1999) 170. 

[34] A.H. Hoang and T. Teubner, Phys. Rev. D60 (1999) 114027. 

[35] A.A. Pivovarov, Phys. Lett. B475(2000) 135. 

[36] A.A. Penin and A.A. Pivovarov, Nucl. Phys. B550 (1999) 375. 

13 



[37] K.Melnikov and A. Yelkhovsky, Phys. Rev. D59 (1999) 114009. 

[38] A. Pineda and F.J. Yndurain, Phys.Rev. D58 (1998) 094022. 

[39] M. Voloshin, Int. J. Mod. Phys. A10 (1995) 2865. 

[40] M. Beneke and V.M. Braun, Nucl. Phys. B426 (1994) 301; I.I. Bigi et al, Phys.Rev. D50 
(1994) 2234; A.H. Hoang, Phys. Rev. D59 (1999) 014039; K.G. Chetyrkin, M. Steinhauser, 
Phys. Rev. Lett. 83 (1999) 4001; S. Groote and O. Yakovlev, |hep-ph/00081~56 . 



[41] Particle Data Group, Eur. Phys. J. C3 (1998) 1. 

[42] V.S. Fadin and V.A. Khoze, Pis'ma Zh. Eksp. Teor. Fiz. 46 (1987) 417; Yad. Fiz. 48 
(1988) 487. 

[43] V.S. Fadin and V.A. Khoze, Yad.Fiz. 93 (1991) 1118; I.I. Bigi, V.S. Fadin and V.A. Khoze, 
Nucl. Phys. B377 (1992) 461; I.I. Bigi, F. Gabbiani and V.A. Khoze, Nucl. Phys. B406 
(1993) 3; J.H. Kuhn, E. Mirkes and J. Steegborn, Z. Phys. C57 (1993) 615. 

[44] B. Kamal, Z. Merebashvili and A.P Contogouris, Phys.Rev. D51 (1995) 4808. 

[45] A.A. Pivovarov, Phys.Rev. D47 (1993) 5183. 

[46] K.A. Ispiryan et al, Yad.Fiz. 11 (1970) 1278. 

[47] F.J. Yndurain, [hep-ph/00073331 , [hep-ph/0008007 . 



[48] M. Jezabek, J.H.Kiihn, M.Peter, Y.Sumino and T. Teubner, Phys. Rev. D58 (1998) 
014006. 



14 



